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A method for computing renormalization constants in the Rome Southampton scheme with vol- 
ume sources and arbitrary momenta is described. This new method is found to enable controlled 
and precise continuum extrapolations and opens the way to compute the running of operators non- 
perturbatively in the Rome Southampton scheme. We describe this in detail and exhibit several 
examples of lattice step scaling functions. 
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1. Introduction 

Regularization invariant (RI/MOM) schemes [|T|] are extremely useful for renormalizing matrix 
elements in lattice gauge theory. They are simple to implement for arbitrary operators and there 
are perturbative calculations, [||] [§J, up to three loops relating RI/MOM to MS. Using momentum 
sources, [|J], excellent statistical precision is also obtained. The leading uncertainties are systematic 
effects: spontaneous chiral symmetry breaking and 0(4) breaking. Chiral symmetry breaking leads 
to a p- pole, violation of Ward identities and a strong mass dependence in the vertex functions. 
Subtracting the condensate term is possible [0], however a change in kinematics is sufficient to 
suppress almost all of the chiral symmetry breaking effects [Q]. In [Q] the incoming and outgoing 
momenta are given the same magnitude and different directions so that the momentum transfer is 
non-zero and equal in magnitude to the external momenta. These non-exceptional kinematics are 
sufficient to suppress spontaneous chiral symmetry breaking by a factor \. 

With these improvements 0(4) symmetry breaking is now the main source of uncertainty. We 
propose to use twisted boundary conditions to remove it. This has immediate impact in reducing 
the error on many phenomenologically important parameters where renormalization is a large and 
often dominant contribution to the total uncertanity. In this work we discuss the quark mass [f7|] and 
kaon bag parameter, Bk [Q, §, §], but the technique also has relevance for many RBC-UKQCD cal- 



culations, for example K ^ tztz matrix elements [10], distribution amplitudes [11, 12] and nuclear 



form factors [13]. 

The largest remaining systematic error is perturbative. Usually the simulated momenta must 
satisfy the condition, 

Aqcd«P 2 «(^) 2 - (1-1) 

The lower bound is for convergence of perturbation theory while the upper is in order to have small 
discretization artefacts. We propose to control the discretization error by continuum extrapolation, 
convergence of perturbation theory can then be improved by choosing a high momentum scale 



before matching. To increase the energy scale while avoiding large lattices we will step scale [14] 
and we outline how to do this using RI/MOM renormalization. 



2. Twisted Boundary Conditions 

On the lattice 0(4) symmetry is explicitly broken to H(4). This means the same quantity 
computed using inequivalent momentum directions will have a different Symanzik expansion de- 
pending on the direction. Matching the same physical momentum and the same direction on two 
different lattices with this constraint is generally difficult. In order to chose momenta with arbi- 
trary magnitude in a given direction we use twisted boundary conditions JT^]. Consider Green's 
functions of bilinear operators q(x)rq(x), T is a Dirac matrix. With twisted boundary conditions 
the quark field satisfies, 

q{x + L) =e id q{x). (2.1) 
Let q(x) = e' Bx q(x) with B = This modifies the Dirac operator, 

D = {0+M) ^D = (0+i0+M). (2.2) 
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The inverses are related by S(x,y) = e lB ( x y 'S(x,y). The momentum source method gives us the 
Fourier transformed propagator, G(z,p), where p is a Fourier mode and 

^D(y,z)G(z,p)=e i Py. (2.3) 

z 

G(z,p) is related to G(z,p) , the Fourier transformed propagator, by 

G(z,p) = £e-'^-%(z,xy> v = e- iBz G(z,p + B). (2.4) 

.V 

Thus, solving the Dirac equation with a plane wave source and twist B gives the momentum space 
propagator G(z,p + B). B is arbitrary and so allows arbitrary momentum in any direction. 
To compute vertex functions Ap we use, 

G r (/;i+£i,/?2 + £2)=l75^ I+B,) *^ (2-5) 

*,y 

Il r (p) = (G- l (p l ,p l )G r (p u p 2 )Y 5 [G- l (p 2 ,p 2 )}^5) (2.6) 

A r =^Tr(U r P r ). (2.7) 

Pr is a projector^]. We pick p\ = (— 1 , 0, 1 , 0) p 2 = (0, 1 , 1 , 0) to minimize £ ; pj. 

Figure[I] (a) shows the effect of the twisted boundary conditions on the data. The 0(4) breaking 
scatter is removed. We simulate momenta with one orientation and continuously varied magnitude 
so that the vertex function is a smooth function of p 2 . The fact that we have a valid Symanzik 
expansion then enables us to take the continuum limit in an unambiguous way by choosing the 
same direction on every lattice. 



3. Step Scaling with RI/MOM 

Twisted boundary conditions mean that we can now study the same physical quantity, at the 
same scale with the same Symanzik expansion on any lattice. This enables a controlled contin- 
uum extrapolation of renormalized quantities Zo(0). Using continuum extrapolation to control 
discretization effects the majority of the remaining uncertainity comes from perturbation theory. 

Rather than try to satisfy the window condition [O] on a single lattice we consider a series of 
lattices of successively descreasing volume and use these to step up to high energy. On each lattice 
we still need to compute vertex functions with small discretization error. The lower limit on the 
momenta is now given by the requirement that we do not resolve the finite volume of the lattice. 
This gives the step-scaling window condition, 

4) 2 «v«n 2 (3-D 



aL V a ) 

By simulating a series of lattices with overlapping scaling windows we will be able to compute 
continuum limit step scaling functions. 

Explicitly, on each lattice at a given quark mass we can compute the ratio, 

R {p,a,m) = — r = ^r7 r- (3.2) 

Ao[p,a,m) Z A (p,a,m) 
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(a) Twisted -vs- Untwisted 



(b) Ay and A& 



(c) As and Ap 



(d) A r 



Figure 1: (a) The axial(red) and vector (green) verticies computed at exceptional momentum [|l]] with a 
volume source at fixed quark mass m q = 0.03 on a 16 3 x 32 x 16, /3 = 2.13 lattice, using 10 configurations 
[p^|. The red and green points use Fourier modes while the blue data uses twisted boundary conditions to 
vary the momentum. Twisting removes the 0(4) breaking scatter. (b,c,d) The bilinear vertex functions in 
the chiral limit computed on 16 3 and 24 3 lattices. All parameters except volume are identical. The volume 
dependence is strongest for the axial and vector currents although it is still only a fraction of a percent 



here is the axial vector vertex function. We use it to divide out a factor of the field renormaliza- 
tion Z q in favour of which can be computed more accurately. The chiral limit is 

Zo(a,p) = limZA(p,a,m)Ro(p,a,m) (3.3) 

m-s-0 

All of the scale dependence is in the renormalization constant and so the factor needed to 
change the scale from p to sp, where s is a constant, is 

v / x ,. Rff(sp,a,m) Z {a,sp) 

Z (p,sp,a) = lim -— = — -. (3.4) 

m-K) Kff[p,a,m) Zo{a,p) 

Continuum extrapolating this gives 

Z (sp) ( /"M y(x) \ 

<*o(p,sp) = UmLoip,sp) = 7 , i =exp[ / -nT^dx . (3.5) 

Z [p) \Ja(p) P{X) J 
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A series of n lattices to each raise the scale by a factor s gives 

(O^Gu)) = (0 SMOM Cp)} x o l {p,sp) ... x a^(s"- l p,s n p) 

\ + c SM0M ^ m cc s (n= s n p)\. (3 - 6) 

Hopefully, perturbation theory is well convergent at the last step where the scale is large. 

3.1 Setting the Scale 

In order to set the lattice scale accurately a different definition of lattice spacing for each 
volume is needed. This means that between two successive steps the same momentum s n p will be 
defined in different ways. The two definitions must be such that they agree in the continuum limit 
so that products like Oq{s u ~ 1 p,s n p)Oq +1 (s n p : s" +l p) are well defined. We need a family of scale 
setting quantities {qi(a)} that depend on shorter distances as we reduce the volume between steps. 

n n / \ cont / \ cont / / \ \ cont 

sp _ s n p (q„-\\ ( q n -\\ -ij m ( 1n-\{a) \ ~ « 



q n (a) q c ™{ \ q n J \ q n J a->0 \ q n (a 

we ensure that scales set on lattice n using q n agree in the continuum with scales set using q n -\. 
We consider a sequence of scales, of the same class as the Sommer scale 



r 2 n F{r n )=C n . (3.8) 

The Sommer scale ro takes Co = 1.65. Thus a step scaling scheme with scale factor s can then be 
denned choosing p n = s n p and r n = ^ as follows: 

• Determine a(po,pi) in continuum limit holding ropa fixed such that r^F{rQ) = Co 

r 2 

• Determine C\ = jiF(-j) in continuum limit holding ro fixed 

• Decrease L by ~ - without fine tuning 

• Determine a(pi,p2) in continuum limit holding r\p\ fixed such that r\F{r\) = C\ 

r 2 

• Determine C2 = -kF\-j) in continuum limit holding r\ fixed 

• Decrease L by ~ - without fine tuning 



s 

etc... 



The guideline r < | should ensure finite volume safety and using the tree level improved po- 
tential [ |l9| ] helps reduce discretization effects. Scale setting in this way will be difficult at short 
distances where the potential runs logarithmically however several steps should be possible before 
this. To investigate the volume dependence of the vertex functions themselve we performed mea- 
surements of identical operators on j8 = 2.13 lattices with 16 3 x 32 x 16 [20] and 24 3 x 64 x 16 



[^] volumes. The results are shown in figure [T] (b,c,d). The volume dependence exists but is much 
smaller than the error from scale setting. Although extrapolations to infinite volume are possible 
they are probably not necessary until the error from setting the scale is reduced significantly. 
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4. Results 



Full details of our ensembles have been reported elsewhere [16], here we show step scaling 
functions for quark mass, quark field and Bk and comment on several interesting features. Figure 
^| (a) shows the step scaling function for the quark mass in the SMOM scheme [||]. The perturba- 
tive and non perturbative results are in quite good agreement. Figure || (b) shows the quark field 
renormalization. At non-zero lattice spacing the perturbative running is in the opposite direction 
to the measured running. This highlights the dangers of applying perturbation theory at fixed lat- 
tice spacing. Only in the continuum limit is it guaranteed that perturbation theory will accurately 
describe the running at high energy. 

Figure || shows the step scaling function for Bk in two schemes p7[]. The choice of in- 
termediate renormalization scheme greatly affects the agreement with perturbation theory here 
SMOM — (g, fi) is apparently optimal. Figure ||(b) shows that on our coarser lattice the perturbative 
and non-perturbative running agree better than in the continuum. Again this emphasises the danger 
of having entangled discretization and perturbative errors. 






: 


. 1^(8,3') p=2.23 

• «8,<») 

LO 

— NLO 

— NNLO 


1 i I 1 1 1 1 1 1 1 1 1 1 1 



(a) a z SMOM 



(b) (T Z SMOM 



Figure 2: (a) Step scaling function for Z m in the continuum limit. Details of the continuum limit are in 
[p"6||. In physical units this corresponds to between 2 and 3 GeV. High order perturbation theory describes 
the running quite well, (b) The step scaling function for Z q from 2 to 3 GeV. Note that the continuum limit 
flips the direction of the running to agree, at least in sign, with perturbation theory. 



5. Conclusions 

Volume source, non-exceptional, twisted boundary condition renormalization is seen to give 
MOM scheme renormalization constants with very small statistical and systematic errors. There- 
fore this will be the preferred method of renormalizing matrix elements in future RBC-UKQCD 
measurements. Further, thanks to the excellent precision available, step scaling arbitrary operators 
now becomes feasible. We have shown here some promising preliminary results and are presently 
computing step scaling functions more accurately and for more operators. We believe this can 
greatly reduce the uncertainty in matching lattice calculations to perturbation theory. 
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(a) Z BK , non-exceptional, SMOM—(^,ff) (b) Z BK , non-exceptional, SMOM - (y^i,/^) 



Figure 3: Step scaling function Bk from 2 to 3 GeV in two schemes compared to the perturbative running. 

References 

[1] G. Martinelli, C. Pittori, C. T. Sachrajda, M. Testa and A. Vladikas, Nucl. Phys. B 445 (1995) 81 
[2] K. G. Chetyrkin and A. Retey, Nucl. Phys. B 583, 3 (2000) 

[3] C. Sturm, Y. Aoki,N. H. Christ,T. Izubuchi,C. T. C. Sachrajda,A. Soni, Phys. Rev. D 80 (2009) 014501 
[4] M. Gockeler et al., Nucl. Phys. B 544 (1999) 699 
[5] J. Noaki et al, Phys. Rev. D 81 (2010) 034502 
[6] Y. Aoki et al., Phys. Rev. D 78 (2008) 054510 

[7] D. J. Antonio et al. [RBC-UKQCD], Phys. Rev. Lett. 100, 032001 (2008) 

[8] C. Allton et al. [RBC-UKQCD Collaboration], Phys. Rev. D 78 (2008) 1 14509 

[9] C. Kelly, P. A. Boyle and C. T. Sachrajda, PoS LAT2009, 087 (2009) 

[10] T. Blum et al. [RBC Collaboration], Phys. Rev. D 68, 1 14506 (2003) 

[11] M. A. Donnellan et al., PoS LAT2007 (2007) 369 

[12] P. A. Boyle, M. A. Donnellan, J. M. Flynn, A. Juttner, J. Noaki, C. T. Sachrajda and R. J. Tweedie 
[UKQCD Collaboration], Phys. Lett. B 641 (2006) 67 

[13] Y. Aoki et al. [RBC-UKQCD Collaboration], Phys. Rev. D 78 (2008) 054505 

[14] M. Luscher, P. Weisz and U. Wolff, Nucl. Phys. B 359, 221 (1991). 

[15] J. M. Flynn, A. Juttner and C. T. Sachrajda [UKQCD Collaboration], Phys. Lett. B 632 (2006) 313 

[16] R. Arthur and P. A. Boyle [RBC-UKQCD Collaboration], arXiv: 1006.0422 [hep-lafj. 

[17] RBC-UKQCD "Continuum limit of B K from 2+1 Flavor Domain Wall QCD", in preparation. 

[18] S. Necco and R. Sommer, Nucl. Phys. B 622 (2002) 328 

[19] G. S. Bali and P. Boyle, arXiv:hep-lat/0210033. 

[20] C. Allton et al. [RBC and UKQCD Collaborations], Phys. Rev. D 76 (2007) 014504 



7 



